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We report on a new topological vortex solution in U(l)xU(l) Maxwell-Chern-Simons theory. 
The existence of the vortex is envisaged by analytical means, and a numerical solution is obtained 
by integrating the eqnations of motion. These vortices have a long-range force because one of the 
U(l)s remains unbroken in the infrared, which is guarded by the Coleman-Hill theorem. The sum 
of the winding numbers of an ensemble of vortices has to vanish; otherwise the system would have a 
logarithmically divergent energy. In turn, these vortices exhibit classical conhnement. We investigate 
the rich parameter space of the solutions, and show that one recovers the Abrikosov-Nielsen-Olesen, 
U(l) Maxwell-Chern-Simons, U(l) pure Chern-Simons and global vortices as various limiting cases. 
Unlike these limiting cases, the higher winding solutions of our vortices carry non-integer charges 
under the broken U(l). This is the first vortex solution exhibiting such behavior. 


I. INTRODUCTION 

Vortices are topological defects that were first dis¬ 
cussed in the context of type-II superconductors by 
Abrikosov [T], where the core of a vortex is in the nor¬ 
mal fluid phase whereas outside the core is in the super¬ 
fluid phase. The relativistic generalization of vortices was 
given by Nielsen and Olesen [2] for the Abelian Higgs 
model. Vortices arise in field theories with degenerate 
vacuum manifolds, whose first homotopy group is non¬ 
trivial, 7ri[AI] 7 ^ /. According to Kibble’s classification 
[3], e.g., the degenerate vacuum of a spontaneously bro¬ 
ken U(l) theory has 7ri[U(l)] = Z (see e.g., Refs. [H [S] 
for reviews). 

The Abrikosov-Nielsen-Olesen (ANO) vortex has no 
electric charge, but has quantized magnetic flux, <I>b = 
27rn/e, where e is the gauge coupling constant and n S Z 
is the winding number of the Higgs field corresponding 
to different topological sectors classified by 7ri[U(l)]. As 
both the gauge and scalar fields are short range, they do 
not exhibit long range interactions. 

Interesting vortex solutions accompany the addition of 
a Chern-Simons term mm I IL which 

breaks the P and T invariance of the theory and gives 
a mass to the photon. It was shown in Ref. [5] that 
if a Chern-Simons term is added to the Abelian Higgs 
model, the vortices carry both a quantized magnetic flux 
= 27rn/e and charge Q = where ^ is the Chern- 
Simons coefficient. Similar to the ANO vortex, the inter¬ 
action is short range and the charge is screened as the 
gauge field is higgsed (see, e.g.. Refs. mm for a review 
of various applications of Chern-Simons vortices). 

Generally, Chern-Simons terms will appear in the con¬ 
text of finite temperature four-dimensional gauge theo¬ 
ries such as the standard electroweak theory m- Upon 
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dimensionally reducing from four to three dimensions and 
integrating out the fermions, non-zero Matsubara modes 
of the gauge bosons and the zero Matsubara mode of 
the temporal component of the gauge fields one obtains 
Chern-Simons terms [T2]. They are also used as effec¬ 
tive field theory models to study the quantum Hall ef¬ 
fect [laiii]. Here, we specifically consider U(l )2 x U(l)^ 
theory with a Chern-Simons mixing term, as given by 
the action Q. In fact, the Chern-Simons mixing term, 
in Q is the 2 -b I dimensional version of 
the BE theory m- 

In this work, we report on a new class of vortex so¬ 
lutions in U(l )2 X U(l)^ Maxwell-Chern-Simons theory. 
One of U(l)s is spontaneously broken by a complex scalar 
field, whereas the other remains unbroken. As a result, 
the new vortex is charged under the unbroken U(l)^, in 
addition of being charged under the broken 0(1)^, and 
it mediates a long-range force. Therefore, an ensemble 
of vortices and antivortices will be confined to minimize 
the energy of the system. This is the dynamical realiza¬ 
tion of the classical conhnement that was pointed out 
by Cornalba and Wilczek m and de Wild Propitius 
m Since our vortices carry magnetic huxes, they will 
also exhibit non-trivial statistics in the infrared. Thus, a 
collection of these vortices will behave like anyons with 
long-range helds. We also show that the parameter space 
of these vortices is vast and includes the limiting cases 
of various known vortex solutions: ANO, U(l) Maxwell- 
Chern-Simons, U(l) pure Chern-Simons, and global vor¬ 
tices [T5]. Interestingly enough, we Hnd that unlike these 
limiting cases, the U(I)j, charge and the U(l)^ magnetic 
hux of the higher winding solutions of our vortices are not 
integers times the charge and flux of the lowest winding 
solution. This is the first vortex solution exhibiting this 
behavior. 

It is crucial that the model we consider does not have 
a self Chern-Simons term, for the U(l)^ 

gauge field which would otherwise spoil its long-range 
behavior. Then, one wonders if quantum corrections can 
generate such a term that destroys the nice long-range 
property of the vortices. Fortunately enough, if this term 
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is absent on the tree and one-loop level, which is the case 
at hand, then the Coleman-Hill theorem [19] guarantees 
that this term will not be generated at any higher loop 
level (see also Ref. El). 

The plan of this paper is as follows. In Sec. [IT] we in¬ 
troduce the Chern-Simons theory with the Chern-Simons 
mixing term for both the mixed and unmixed basis, and 
then discuss their basic properties. In Sec. m we give 
a proof of existence for topological vorticies that are 
charged under the long-range U(l)^, and then present our 
results for the numerical solutions. In Sec. |IV[ we calcu¬ 
late the flux, charge and energy of the vortex solution and 
of a vortex-antivortex pair. We then show that the en¬ 
ergy of the vortex-antivortex system is finite whereas the 
single vortex energy is logarithmically divergent, hence 
the vortices are classically confined. We conclude with a 
summary of our results and discussion in Sec. |V| 


II. TOPOLOGICALLY MASSIVE U(l)xU(l) 
THEORY 


We consider two topologically massive Abelian gauge 
fields 3^^ and with corresponding gauge groups 
U(I)^ X U(1)y, and a complex scalar field (/? that is cou¬ 
pled to a linear combination of and 


S = 




OL 
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- + |(a^ - igiy^ - mWM (1) 


theory as follows: e = \/gi -f gi = 2/iw tan 0, and 
fi 2 = 2/rw(tan^ 6 — 1). The Chern-Simons coefficients gi 
and /i 2 as well as the parameter A have mass dimension 
M, while the coupling constant e and the vacuum expec¬ 
tation value V have mass dimension We set c = 1, 

h= 1, = 1, and use the metric = diag(l, —I, —1) 

in what follows. 

It is a simple exercise to study the fluctuations about 
the vacuum \{p\ = u in Eq. Q. First, the gauge field 
Api carries a single massless degree of freedom, thanks 
to the unbroken U(l)^. Writing the complex field g) as 

= {v + we find that there is a single massive 

radial field h in the infrared. In addition, the would- 
be Goldstone boson, 11, is eaten by the massive 
field. In fact, the mass of receives contributions from 
three sources: the self Chern-Simons term /i 2 , the Chern- 
Simons mixing term gi, and the Higgs vacuum expecta¬ 
tion value. This will be clear from our vortex solution, 
as is evident from Eq. 0 below. Thus, the field has 
two degrees of freedom, and we recover the total sum of 
the four degrees of freedom we started with. 

One wonders whether the condition ([^ and hence the 
spectrum described above, especially the massless U(l)^ 
field, are not spoiled by quantum effects. In fact, a one- 
loop calculation in the theory described by Q does not 
yield any corrections to the Chern-Simons terms m- Be¬ 
sides, according to the Coleman-Hill theorem m, there 
are no more corrections to these topological terms other 
than the one-loop contribution. Therefore, the massless 
U(l)^ gauge field is protected against quantum effects. 


III. CHARGED VORTEX SOLUTION 


where = dfj.y„-d^yf_, and For 

generic values of and , each of the fields y^ and 
has a single degree of freedom which is screened in the 
infrared, thanks to the topological masses. Adding the 
two degrees of freedom of the complex scalar, our system 
has four degrees of freedom in total. In the following it 
will be useful to go to the new basis and Z^: 

y^ = cos OA^ + sin 6Z^ , 
yy^ = - sin OA^J, + cos , 

where tan9 = gxjg^- Now, we fix 

~~ tan 9 (1^) 


Throughout this work, we seek cylindrically symmetric 
vortex solutions of the theory given by the action ©.By 
varying the action [Eq. Q], we obtain the field equations 

+ = 0 , 

dfiZd'^ + + g2t^"^Zp^ + j" = 0 , ( 5 ) 

DpD^ip + ^ -v^)(p = Q, 

where we defined the current as 

f = . ( 6 ) 


to obtain the action for the corresponding U)!)^ x U(l)^ To this end, we take the cylindrically symmetric Nielsen- 
theory Olesen like Ansatze, namely. 


S = 






( 4 ) 


where ^11 — d^Au duA^-, Z^^ — d^Zjj duZ^^ and 
= d^ — ieZ^. The parameters of the U)!)^ x U(l)^ 
theory are related to the ones in the 11(1).^,^ x 11(1)^ 


= vf{r)A'^^ , Zj = , 

er^ 

Zo = eZo(r), A = (7) 

er^ 

A = eAo(r), 

where the profile functions /(r), Z(r), Zo(r’), A(r), and 
Ao(r) are dimensionless. Using these Ansatze in the equa- 









3 


tions of motion ([^, we obtain 

- (n - + e^Zlf - ^(/2 _ 1)/ = 0, 


r 

Z’ 


Z" - — + 2e^v^f{n - Z) - 2e^r{^JLlA'Q + /^2^o) = « > 

Z" + ^ - 2e^v^fZo - + f,2Z') = 0, 

/V, C^T 


12 

r 

A' 


( 8 ) 


A" - 2/iie^rZ' =0, 

T 

A'' + ^-‘^-^Z' = Q. 


Note that in the limit /xi = 0 the two U(l) sectors decou¬ 
ple and we obtain the equation of motion for the normal 
11 ( 1)2 Chern-Simons vortex [8]. The last two equations 
in Eq. ® can be integrated to find 


A — 2^\(^tZq T^it ^ 


A' =^Z- 
e^r 


V. 


(9) 


where I?i and I ?2 are integration constants. In order to 
determine the constants I?i and I ?2 j we examine the near 
core and large r behavior of the system. As we shall show 
the behavior of Z(r) near the core goes like 


III A 


in Sec. 

r^, and hence, one has to set I ?2 = 0 in order to have 
a regular solution of the electric field eAQ(r) at r = 0. 
Besides, a regular solution for A'{r) at large r demands 
that Vi = —2/iie^Zo(oo). However, since nonzero Zo(oo) 
leads to a quadratically divergent energy [see Eq. (27)], 
it has to be set to zero, so does Vi. 

Before delving into the detailed vortex solution, one 
can read the physics of the vortex solution from the 
second equation in (1^. This relation states that start¬ 
ing with a single U(l )2 Chern-Simons vortex, i.e. setting 
/Ti = 0, which has an asymptotic Z solution of the form 
Z{oo) = n ,n G Z, and turning on a small will cause 
the vortex to acquire a long-range electric field propor¬ 
tional to fill 


Ea = eAp ^ 


2/Xin 


er 


( 10 ) 


Therefore, our vortices will carry a long-range field, 
thanks to the unbroken U(I)^. This physics will be con¬ 
firmed by detailed analytical as well as numerical checks, 
as we show below. 


Eq. (IS, setting / = I, and neglecting the derivative and 
©(l/p) terms: 


c. u u 

- g2^2 + 2fj.l ’ 


Zo{oo)=0, /(oo) = I, 

( 11 ) 


A(oo) = nC 2 , Ao(oo) = 


2fiiv‘‘n 

+ 2^1 Cl 
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, e^r 
in — 


where the constants Ci ,2 are determined numerically. Un¬ 
like the single U(I )2 Chern-Simons vortex, the Z wind¬ 
ing is not an integer. This peculiar behavior will have 
its dramatic consequences on the single as well as the 
multi-vortex solutions as we discuss later on. In the limit 
0, we find Z{oo) = n as expected for the single 
U(I)z Chern-Simons vortex. 

Near the core, r —>• 0, the fields can be expanded 
in a Taylor series with arbitrary parameters. Requir¬ 
ing that the physical fields are continuous at the ori¬ 
gin and forcing the expansion to fulfill the second order 
field equations ([^ , we can fix all but the five parameters 
OOO) -ZOO) Z2,a2, fl- 


fir) = + 0(r3), 

Z{r) = ez 2 r^ -b O(r^), 

+ ,j2) 

e 

Air) = ea2r^ + C>(r"‘), 

Aair) = -h Oir ). 

e 


Note that the value of ago is purely a gauge choice, which 
doesn’t contribute to the equations of motion (§. Below 
we will arbitrarily set ago = 0. 


B. Existence of the Solution 


A. Boundary Conditions 

In the usual Nielsen-Olesen vortex solution, the asymp¬ 
totic behaviors of the fields are determined easily by their 
regularity at the core and the finiteness of the energy; 
namely, the profile functions vanish at the core, Higgs 
goes to its expectation value, whereas the gauge field 
goes to a pure gauge value determined by the vanish¬ 
ing of the kinetic energy of the Higgs field at infinity. In 
our model, the boundary conditions of the profile func¬ 
tions at r —>• 00 are obtained by substituting Eq. ([^ into 


Before moving to the numerical solution of the field 
equations, in this section we sketch a proof of existence of 
the vortex solution. To this end, we analyze the system of 
equations ([^ at asymptotic infinity, r —> 00 , taking into 
account the first order equations ([^. Since our vortices 
carry a long-range U(I) field, it is expected that the far- 
held will follow a power-law behavior. The asymptotic 
power law behavior can be obtained by expanding the 
prohle functions as 'k(r) = 'k(oo) -b J2m=i'f’m/r'^ and 
solving for using Eq. ([^. Upon performing this ex¬ 
pansion, all constants tpm can be hxed and we get for the 
asymptotic prohle functions d>(r) —>■ tE'°°(r) to leading 
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order: 

rir) 

Z°°{r) 


ZS°ir) 


A^{r) 

A°°(r) 




1 - 




+ 0(l/r4), 




^(«^) - U2T'r2\4 2 + ^(lA") , 
Z2n^nl^i2 , / 6 


A(e2r:2 + 2^^)^ r'^ 

16nVI 


^o(oo) + 
nC2 + 


A(e^v^ + 2^^)"^ 

32e^n^ fil^2 


0{l/r^), (13) 

0(1/A), 


A(e2D^ + 2fil)^ 


0 ( 1 //). 


Note that in the ^2 = 0 limit, ^o(^) = Al(r) = 0 to all or¬ 
ders in the large r expansion. In addition, in the —> oo 
limit, all the profile functions except /(r) vanish identi¬ 
cally, hence we recover the U(l) global vortex solution 

m- 

Now we expand the profile functions around the large 
r limit as 


f = r+Sf, Z = Z°^+ SZ, Zo = Z^ + SZo , (14) 


and plug these Ansatze back into the equations for the 
profile functions Q. To the leading order in the fluc¬ 
tuations Sf,SZ,SZo (and neglecting terms of the form 
{6f, 5Zo,5Z)lr‘^), we get 

Sf" + ^ _ Xv^Sf « 0 , 
r 

fs7' 

6Z" -(2e\2 + i^il)6Z - 2e‘^rfi2SZ' « 0 , (15) 

r 

R7' 9 

bZl + ^ - (2^v^ + AnDdZo - -^^125Z' « 0 . 


By first writing 5Z(r) = e^/rzir) and SZq^t) = 
zo{r)/{e^/r), then solving for z{r) and 2 ;o(r) at large r, 
we find following regular solutions to the homogenous 
equations Eq. (15) 


Sf(r) ~ 

ey/r 

SZ{r) ~ Cae/re”^^’’, (16) 

5Zo(r) cs 

eyjT 


where AI^ is the mass: 




= -A 2 I + inl + 2e'^v'^ + . 


(17) 


\og{e'^v^/2fil) 


The solution pT| describes the intermediate region 
> r, where the equality sign results by 


V Xv 


equating / from the power-law behavior (13) to the expo¬ 


nential one (16). Assuming that a vortex solution with ra¬ 


dius Tc exists, then the profile functions and their deriva¬ 
tives have to be continuous at re- Thus, we match the 
profile functions and their derivatives in the small r limit 


given by Eq. (12 1 with the large r limit in Eq. (14). If 


the solution exists, then we should have the right number 
of free parameters. Notice that if we use the first order 
equations 0 instead of the second order ones ([^, we 
obtain a non-trivial relation between 02 and zqq : 


02 — ^J-lZoo , 


(18) 


which reduces the number of the free parameters by one. 
We are then left with three free parameters at the core 
-Z 00 j-Z 2 ,/i, and five free constants at infinity Ci_ 2 , 3 , 4 , 5 - 
Now matching / and its derivative across Tc gives /i and 
C 5 , Z and its derivative across Vc gives Z 2 and C 3 , and Zq 
and its derivative across Tc gives zqo and C 4 . Matching 
Ag and A! across does not give new information since 
both of these functions are dependent on Z and Zq, as 
is clear from Eq. Finally, we can solve for Ci _2 by 
matching A and Aq across Vc- 

The explicit expressions of Ci 3 4^5 are cumbersome and 
not very illuminating, and we refrain from giving them 
here. It will turn out that the value of C 2 determines the 
magnetic flux of and partially the electric charge of 


Now two comments are in order. First, one can read 
the masses of the particle spectrum from Eq. (16): the 
field is massless, the Higgs mass is -x/Au, while the 
mass is given by Eq. ( [T7| ). The mass gets contribution 
from the Higgs vacuum expectation value, after eating 
the would-be Goldstone boson, and from the topologi¬ 
cal Chern-Simons terms. Thus, as stated before, the Z^ 
field has two degrees of freedom. Second, we note that 
in the limit /i 2 = 0, the coefficients zqo, C 2 and C 4 are 
identically zero for all values of /ii as we checked numeri¬ 
cally. In this limit, the Chern-Simons vortex degenerates 
to Abrikosov-Nielsen Olesen vortex. The vanishing of C 2 
for /42 = 0 means the absence of the A^ magnetic flux 
and Z^ electric charge as we detail below. 


C. Numerical Results 


The full numerical solution of the second order equa¬ 
tions (|^ is obtained via a shooting method as imple¬ 
mented in |20j . Starting from a small but non-vanishing 
radius r„iin = 10“^, in units of e^, and using the small r 


expansion ( 12 ) to sixth order, the shooting method finds 
the value of the four free parameters that lead to a solu¬ 
tion satisfying the first four boundary conditions at large 
distance rmax as given in Eq. 0 - Here, to reach the 
boundary conditions to an absolute precision of 10 “® at 
?'max = 27, we use 128 digits of precision to solve the set 
of non-linear differential equations (|^ . Note that the last 
boundary condition at infinity in Eq. (11) is satisfied au¬ 


tomatically. The profile functions, electric and magnetic 
fields, charges and energy density for the specific case 
?; = 1, A = 1, = ^2 = 1/4 are shown in Figs. 

and for different winding numbers, n, as a function of 
r, which is in units of e^. 

The existence of a well defined solution is ensured by 
checking that it follows the expansion at small r given in 
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FIG. 1. From top to bottom, the profile functions for n = 1, 
n = 2, n = 3 vortex solutions as a function of r in units of e^. 
Far from the core, all the fields converge as power law to a 
constant value except for Ao which increases logarithmically. 
At large radius, the functions /, A, Z, Zq, and Ao tend to 
their asymptotic values as given by Eq. fig. 


FIG. 2. The electric and magnetic fields as given by Eq. (191 
for n = 1, 2, 3 from top to bottom. All fields decay as power 
law outside the core except for Ea, which is a long-range field. 
The thick green curve is the fit to the asymptotic electric field 
given by QA/r. 


Eq. (12) and merges smoothly to the asymptotic behavior 
given in Eq. (13). In Fig. we show how the profile 
function /(r) converges towards 1 at large distance and 
see that it satisfies the large r expansion [Eq. ([T^]. As a 
further test we checked that the solution does not depend 
on the value of and Cmax as long as they remain small 
and large enough, respectively. 


tions using Eq. ([^: 

Ez=eZ'o, Bz = = — , 

2 er 

Ea = eA' , Ba = = - , (19) 

2 er 

Zo = eZo, \D^\^=v^f^ + v\n-Zf^-^. 


IV. PHYSICAL PROPERTIES OF CHARGED 
VORTICES 


A. Magnetic Fluxes and Charges 

The electric and magnetic fields as well as the kinetic 
term for the Higgs held can be related to the prohle func- 


The magnetic hux of and helds for the charged 
vortex solution are 


^Bz = 


^Ba = 



Z-Al = 



Z{r) 

er 


Stt 

rdO = — Z{oo). 


Adi = 



Mr) 

er 


Stt 

rdO = —A(oo), 


( 20 ) 


where is a circle enclosing the vortex at inhnity. Using 





























6 




FIG. 3. Zooming in for the large r behavior of the Higgs profile 
/ to monitor the convergence of 1 —/(r) (red full line) towards 
zero. At distances larger than the vortex core, the fields first 
decay exponentially (the green dashed line) as in Eq. (16l. 


At larger r, one reaches the small power law tail (blue dotted 
line) of order /rf as given in Eq. (131. Note that the constant 
C 2 is fitted b ut t he constant ci is known analytically from 
equation Eq. (131. 


Eq. (11), we find 


= 


27rn 


2 2 
e V 


e + 2 /xf 


= ( 21 ) 


Although we have used the asymptotic values of the fields 
to calculate the fluxes, it should be clear that these fluxes 
originate from the near-core region of the vortex since 
both and A/j, are screened outside the core, as is clear 
from Eqs. (16 1 and (13) and Fig. As we discussed be¬ 
fore, C 2 and hence vanish identically at /i 2 = 0 for 
all values of fii. Note that in the /ri —>■ 0 limit Eq. (21) 
reduces Xo ^Bz = ‘i-Tmje as expected for a single UfTJ^ 
Chern-Simons vortex [5]. Unlike the single U(l )2 Chern- 
Simons vortex, the flux in our case is not an integer times 
27r/e, which is attributed to the mismatch between the 
Z,, and 09 windings. More on this point will be discussed 
in Sec. HV^ 

The charge of the vortex under the held can be 
obtained from Eq. ^ as (this is the Noether’s charge) 


Qz = 


/ P‘2'K nOO 

= 2e^u^ / da dr r/^Zo. (22) 

Jo Jo 


Using the third equation in (1^, we have 


Qz = 27re / dr 
Jo 


(rZ')'--(/iiA' + M 2^0 


= 27re 




(23) 


and then, using the boundary conditions in Eq. ^,Qz 
reduces to 


47rn 

Qz =- 

e 


/X1C2 -f /J .2 


-b 2/rj 


(24) 


FIG. 4. The charge density under U(l )2 as given by the inte¬ 
grand in Eq. (22 I for n = 1, 2, 3. Notice that the most of the 


contribution to the charge comes from the near-core region. 


Remembering that C 2 = 0 at /i 2 = 0, we see right away 
that the Noether charge of the vortex under the Z^ held 
vanishes in this limit. In fact, the absence of the Z^ elec¬ 
tric charge and magnetic flux at /r 2 = 0 is not a 
coincidence. We can understand this observation as fol¬ 
lows. For fJ^i = fJ .2 = 0, we recover the normal Abrikosov- 
Nielsen-Olesen vortex which carries only Z^ magnetic 
flux. Turning on a non-zero value for ^ 2 ) keeping /ii = 0, 
the Chern-Simons term will induce an electric charge for 


Z^ held. Recalling Eq. (|22|) — which determines the 


the 

electric charge as a function of Zq — and the discussion 
after Eq. (16), the values of the coefficients zqq and C 4 


are determined upon matching the near-core and the far- 
region values of the profile function Zq and its derivative 
across the vortex wall r^. For /i 2 yf 0 both C 2 and zqo 
are non-vanishing, and hence we obtain non-zero values 
for Qz as easily seen by taking the —>■ 0 limit in Eq. 
(24), which reduces Qz to — 47 rn^ 2 /e as expected for the 
single U(l )2 Chern-Simons vortex. 

Now, let us turn on a non-zero value for fj,i keeping 
^2 7 ^ 0. The first equation in © relates A' to Zq at all 
values of r. At the core, non-zero values of Zq, which 
are expected for ^2 7 ^ 0 , will induce non-zero value for 
the A profile which will induce magnetic field Bj^, and 
hence magnetic flux Thus, we see that the vanish¬ 
ing of ^2 means the vanishing of Zq, and hence Qz, and 
in sequence the vanishing of the flux $■ The interest¬ 
ing feature of Qz is that it is not quantized. Although 
the charge neutrality condition is always satisfied, i.e., 
Qz{-n) = -Qzin), we find that Qzin) ^ nQz{n = 1 ) 
if /ii yf 0 and /i 2 7 ^ 0. If /r 2 = 0 then Qz vanishes, and if 
/ii = 0 the charge is quantized as is obvious from formula 
(24). In fact, the second term of Eq. (24) is directly pro¬ 


portional to n but the first one is not as the asymptotic 
value C 2 depends non-trivially on |n|. This is the first ex¬ 
ample of a vortex with such behavior. The Z charge of 
the vortices with winding n = 1 , 2 ,3 are given in Table |T] 
for ^2 = 1/4, e = V = \ = 1 and several values of fii. 

Far in the infrared, the Z^ field is screened while the 
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n=l 

n=2 

n=3 

/ri = 0 

3.1415 

6.2831 

9.4247 

pi = 0.25 

2.2908 

4.2625 

5.9838 

pi = 0.50 

1.0752 

1.8116 

2.4765 

pi = 1 

0.2091 

0.4066 

0.5672 


TABLE I. Charge Qz for ^2 = 1/4, e = A = w = 1 and 
different values of and winding numbers n. Note that it 
is quantized for fii = 0 (local U(l )2 vortex) and vanishes for 
^ 00 (global vortex), but in general this is not the case. 
The error estimate of the charges is < 0.01%. 


Ao field is long-range, thanks to the unbroken U(l)^. 
The charge of the vortex under U(l)^ can be obtained by 
integrating the electric field over the plane and using 
Stokes’ theorem. Substituting the asymptotic value of Ao 
in Eq. 0 we obtain 



• d£. = 



A'^rde = 


e^v'^ 2^\ 


(25) 


This is the electric charge of the vortex under the U(l)^ 
field as defined from Gauss’s law [note that < 5,4 = 
where given by Eq. (21)]. Therefore, an 
external probe with test charge Qtest will experience a 
forc^ QtestQ^/?'- The charge Qa is zero at fii = 0, in¬ 
creases to a maximum value of 27m/e at /rf = e^r!^/2, 
and then decreases as 2Trnev‘^/fii for 2fil The 

decrease of the electric charge for large values of fii can 
be understood from the equations of motion ([^. For large 
values of fii we can neglect the kinetic terms compared to 
the topological one. Since the kinetic terms are responsi¬ 
ble for mediating the long-range U(l)^ force, we expect 
this force to be suppressed for small values of the kinetic 
terms (see also Fig. [^. 

Before ending this section, let us also note that if the 
charge were defined from Gauss’s law, as we did for 
we would find zero Z^ charge since the Z^ field is 
screened. 


B. Energy of the Vortex 

The energy of the vortex can be calculated starting 
from the Hamiltonian density of the theory defined by 
Eq. ([^. Because the vortex is static, one can alterna¬ 
tively use the Euclidean version of this action. In fact, 
since the Chern-Simons terms do not depend on any 
background metric, these terms do not contribute to the 
energy-momentum tensor and hence to the Hamiltonian. 
In turn, one does not expect that both the Hamiltonian 
and the Euclidean action to have the same functional 
form. Irrespectively, we checked that both formulations 
give the exact same answer for the vortex energy. 
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The interaction between vortices will be elucidated in Sec. 


IV C 



FIG. 5. The energy density of a vortex with n = 1, 2, 3. The 
blue dash-dotted lines show the asymptotic behavior given by 
the derivative of the infrared contribution to the energy given 
by Eq. (291 with respect to r. 


The Hamiltonian density is given b}0 


■Hv 


i {El +BI+EI + B\) + 

+ |D(^P + ^(M2-^2)2. 


(26) 


Next, we split the total Hamiltonian density into three 
parts: the core 'Ha electric T-Lea Goldstone 'Hq. Us¬ 
ing Eq. (IM, 


He 


+2uY^ + ^(f-l)2 , 


nG = v^{n-zr^. 


(27) 


Integrating each term over , the total energy of a vor¬ 
tex is obtained as 


— £c + £ea + • (28) 

The core contribution to the energy, £c can be found nu¬ 
merically (see Fig. [^, while in the IR only the electric 
£ea Goldstone Eq contributions remain, and the en- 


^ Note that our formula differs from the one given in Ref. CQl by a 
sign of the term, which should be positive. Besides, the 

Chern-Simons term that they have in the Hamiltonian density 
should not be included as it does not contribute to the energy 
momentum tensor. 





















ergy can be calculated analytically: 




'Ea 


+ 


^^Z((X))^ + 2v‘^{n — Z{oo)y 


rRiR 


d\nr 


47rn^n^u? , i?iR 
■ In- 


+ 2/i 


(29) 


where we imposed an IR cutoff at i?iR. Exactly like an 
electric charge in 2+1 dimensions, the energy of a single 
vortex is logarithmically divergent. Note that there are 
two contributions to the energy. The first one comes form 
the electric field £ea^ which is expected from the long- 
range U(l)^ outside the vortex core. The second con¬ 
tribution, the Goldstone energy, Sq is more interesting. 
In a trivial background, + = 0, we can study the spec¬ 
trum of the symmetry breaking of a global U(l) symme¬ 
try by writing (p(x) = \v + h{x)] where h is the 

Higgs boson and the phase H is the Goldstone boson. 
This Goldstone boson is a physical massless degree of 
freedom that exists in the spectrum of the theory. How¬ 
ever, once we gauge U(l) this Goldstone will be eaten 
by the corresponding gauge field which in turn acquires 
a mas^ In fact, the physical spectrum does not con¬ 
tain any Goldstone boson since there is no symmetry to 
break: the gauged U(l) is a redundancy rather than a 
genuine symmetry. If the background is non-trivial, as 
in our vortex case, then one has to repeat the same ar¬ 
gument in the given background. In this case, we write 
+(x) = n [/(r) +/i(x)] where /(r) is the 

profile function of the vortex. We also write the field 
as Z^(x) = Z^ix) + Z^{x), where the background solu¬ 
tion Z^ can be read directly from the profile functions 
Z{r) and Zo{r) given above. Now the square of the co¬ 
variant derivative gives [dfj,{n6 + H) — e(Z^ + Z0] . 
Again, one can use an appropriate gauge transformation 
to kill the fluctuating Goldstone field H, thus interpreting 
as the mass of in our non-trivial background. 
What remains is the Goldstone background cont ribution 
[ndf^O — eZ^Y which gives 'Hq in Eq. (271. As we 
showed above, this Goldstone background energy is log¬ 
arithmically divergently Therefore, unlike the Abrikosov- 
Nielsen-Olesen vortices or their single U(l) Ghern-Simons 
cousins, our vortices are not genuine solitons. This forces 
us to consider an ensemble of an equal number of vor- 


® This can be seen once we compute the square of the covariant 
derivative which gives — , where zf^ is the 

fluctuation field. Using the gauge transformation Zfj, —)■ Zfj, + 
dfitl/e (unitary gauge), we immediately recognize 2e^i;^ as the 
Z-mass. 

^ In fact, one can obtain the same behavior in the presence of a 
non-dynamical static magnetic Bz field with non-integer flux. In 
this case, we also find that the Goldstone background energy is 
logarithmically divergent because of the mismatch between the 
winding of the Higgs and gauge fields. We thank T. Sulejmanpa- 
sic for emphasizing this point. 


tices and anti-vortices. These in turn will form vortex- 
antivortex confined pairs that lowers the total energy of 
the system. The next section is devoted to the study of 
the interaction energy of the vortex-antivortex pair. 


C. Interaction between Two Vortices 


In this section, we show that the total energy of a sys¬ 
tem of a vortex and an antivortex pair is finite. Treat¬ 
ing the 2-vortex system is rather complicated. However, 
it will be sufficient to approximate the system as a su¬ 
perposition of a pair of vortices with opposite winding 
numbers located at a large separation R. As long as the 
separation is large enough, the individual solutions do 
not receive a considerable modification by the presence 
of the other vortex. With these assumptions, the total 
scalar and gauge fields of a vortex-antivortex system at 
|x| ^ Tc can be approximated as follows: 


^0 = 0, 


57 ^ 

A ~ *3.4, 

A = 

ZTT |X — X2 


(x-Xi)j (x-X 2 )j 


X - Xi| 
X-Xi| 


X-X2 


(30) 


A 

e 


(X-Xl)j (X-X2) 


|X-Xl| 


|X-X2| 


where, is given by Eq. (25), xi and X2 are the lo¬ 
cations of the vortex and antivortex, respectively. Upon 
substituting Eqs. (30) in Eq. (26), the total Hamiltonian 
density in the IR can be obtained as: 

, SnV^? (R/2)2 

H^n{R,r,a) + r^rj ’ 

where we defined r = |x|, 

ri = yjr'^ — Rr cos a + (77/2)^ , 
r2 = \/r'^ + Rr cos a + (77/2)^ , 


(31) 


(32) 


R is the separation distance between the two vortices, 
and a is the polar angle. Then, the total energy can be 
found upon integrating 'Hi^(r,a) over the polar angle a 
and radial variable r: 


p p‘l'K 

£i^{R) = dr r / da'HiB.{R,r,a) 

Jc Jo 

TT ^ 

+ 2p,l 2 RJr"^ — r^ 


(33) 


Since R ^ r^, the total energy of a vortex-antivortex 
system can be simply written as: 


fia(i?) 


Sirn’^v’^yi R 

+ 2fj,l Tc 


0{rjR). 


(34) 


Here, for simplicity we took the region of integration for 
the radial variable to be C = {0 < r < Rf2—rc and i?/2+ 
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’’c ^ ^ 00 } to remove the contribution of the cores of 

the vortices. It would be more accurate to cut out just 
two discs of radii r<. centered at the cores of the vortex 
and antivortex. However, removing the contribution of 
this tiny strip of radius 2rc only leads to an error of order 
Tc/i?, which is a subleading effect that we ignore here. 

The total energy takes an even simpler form in the 
2^\ <C limit, i.e., when the contribution to Z- 
mass is small compared to the mass coming from the 
spontaneous breaking of the U(l )2 symmetry: 

In^. (35) 

In other words, the interaction of a vortex and an an¬ 
tivortex in this particular limit is exactly like that of two 
point particles with opposite charges in 2 -|- 1 D. 

In the opposite limit 2^\ the total energy 

takes the form 

^iR(^)l 2 /.f»e 2«2 - 47rn\^ In — . (36) 

In this limit, ~ 0, hence, the contribution is mostly 
due to the Goldstone background. Note that this is noth¬ 
ing but the energy of a global vortex-antivortex pair (see 
FigQ. 

To summarize, a system of a vortex-antivortex pair has 
a finite energy and is logarithmically conhned. The total 
energy gets contributions both from the electric held and 
Goldstone background of the vortices. 


V. SUMMARY AND DISCUSSION 


In this work, we have obtained a new vortex solution in 
the U(I) 2 xU(I)^ Chern-Simons gauge theory. These vor¬ 
tices are classihed by a topological number n gZ. Inside 
the core of the vortex the Higgs held is in its symmetric 
phase and both U(l)s are topologically massive. Outside 
the core, the Higgs held gets a vacuum expectation value 
causing the spontaneous breaking of 0 ( 1 ) 2 . The result¬ 
ing Goldstone boson is eaten by the corresponding gauge 
held, namely the Z^ held, which now acquires an extra 
degree of freedom. Thus, the mass of the Z-boson gets 
contributions from both the topological terms and the 
Higgs vacuum expectation value, as is evident from Eq. 
(dzi)- In addition to the massive Z-boson, the vortex me¬ 
diates a long-range force outside its core, thanks to the 
unbroken U(l)^. This adds up correctly to the number of 
degrees of freedom (d.o.f.): inside the core the Higgs has 
2 d.o.f. and each U(l) has a single d.o.f., and outside the 
core the massive Z-boson has 2 d.o.f. while the yl-boson 
(photon) and Higgs each has a single d.o.f. 

Our vortices are charged under the U(l)^ and hence 
two vortices will have a logarithmic interaction due to 
the massless U(l)^ field. The electric field E _4 is given 
by 



r 


47rneu^/ii 
-I- 2/rf ’ 


(37) 



FIG. 6. The field lines of a vortex (left) and an antivortex 
(right). The radial (red) and circular (blue) field lines cor¬ 
respond to t he e lectric field [Eq. (371] and Goldstone back¬ 
ground [Eq. (39 1 ], respectively. 


where is a unit vector in the radial direction. In ad¬ 
dition to the long-range electric field, our vortices also 
interact due to a background Goldstone field. As we 
stressed at the end of Sec. |IVB[ this Goldstone back¬ 
ground is different from the Goldstone fluctuations which 
are completely absorbed by the Z-boson. Let us define 
the Goldstone field: 


G, = vd,e - Z,, (38) 


which can be rewritten as 


G = 



Qg 


-I- 2gL\ 


(39) 


where Bq, is a unit vector in the polar direction, and 
we have defined an effective Goldstone charge Qg (see 
Fig.§. Then, the total interaction energy between vor¬ 
tices with charges {Q\,G}j) and {Q‘\,G'q) is given by 


— (QaQa + QgQg) 111 ~ ■ (40) 


At /Ti = /r 2 = 0 our vortices reduce to the Abrikosov- 
Nielsen-Olesen vortices. Turning on a non zero value of 
fj. 2 , but still setting = 0, we recover the single U(l )2 
Chern-Simons vortex. For both cases the charges Qg and 
Qa are zero, and the vortices do not interact with long- 
range forces. As we turn on a non-zero value for ^ 1 , our 
vortices start interacting logarithmically. For values of 
e^u^/2 the interaction is dominated by the electric 
force, while the Goldstone force is subleading. This pic¬ 
ture is reversed for fj,f e^u^/2 as the Goldstone force 
dominates over the electric one. Both pictures are inde¬ 
pendent of the value of ^ 2 - The only effect of fj .2 is that 
it contributes to the Z mass, as is clear from Eq. 0. 
making it infinitely large as p .2 —>■ 00 . On the other hand, 
taking /ri —)■ 00 , the electric charge as well as the mag¬ 
netic flux <I>b _4 vanish, the later is clear from the fact 
that all the gauge fields decouple in this limit, and the 
interaction is solely due to the Goldstone background. In 
fact, this is the limit where we recover global vortices. Of 
course, in this limit the Z-boson becomes infinitely mas¬ 
sive and decouples. The infinite Z mass can be thought 
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Pure CS vortex 


# 


U{l)zXU(l)^ 

Confined vortices 




11(1)2 ANO Vortex eV /2 




FIG. 7. The parameter space of the vortex solution. The cir¬ 
cle size indicates the core radius, which decreases as we in¬ 
crease both /ii and /i 2 . Our solution covers the whole region 
including the boundaries, and reduces to the known solu¬ 
tions in certain limits of Chern-Simons coefficients fii and fi 2 ■ 
e^u^/2 marks the boundary between electric field dominated 
vs. Goldstone background dominated solntions. 


of as a UV cutoff on the vortices. The parameter space 
diagram of these different limits is illustrated in Fig. 

In Fig. 1^ we show the profile function /, and its asymp¬ 
totic behavior for these regimes. Since the energy of a 
single vortex is logarithmically divergent, its long energy 
tail has to be trimmed by either putting the vortex in a 
container with a finite radius, or by considering an equal 
number of vorticies and antivortices. In the later case, 
the system will lower its total energy by forming con¬ 
fined vortex-antivortex pairs. 

So far, we have not discussed the effect of the mag¬ 
netic fluxes ^Ba ^Bz on the behavior of the vortices. 
It was shown in m that the monodromy of a particle 
{^Bat^Bz) and a remote particle leads to 

an Aharonov-Bohm phase 


exp [f/ri {'^Bz'^'ba + ‘^'bz'^Ba) + ^^J-2<^Bz'^'Bz\ ■ (41) 


This phase gives rise to non-trivial statistics of the vor¬ 
tices, which behave as anyons. One can understand the 
origin of this phase as follows. Given a charge q that cou¬ 
ples to a vector potential An, this particle acquires the 
Aharonov-Bohm phase exp \iq§ = exp[zq$B^], 

where ^Ba lli® Aox of A^^, as the particle makes a 
non-contractible winding around the source Ai^. Now, 
the zeroth component of the equations of motion ([^ (the 



FIG. 8. Profile for the Higgs field 1 — /(r) for different /ri’s 
(M 2 = 0 here). We see that the vortex prohle goes exponen¬ 
tially to the Higgs vev (straight line on the log-linear plot) 
in the case of the ANO vortex {fii = 0) and continuonsly 
changes with increasing fii to a power law, to our vortices 
and to the global vortex in the limit fj,i = 00 . 


Gauss’s laws) can be written in the form 


d xW ■ E_4 = 2 fjLi^Bz 


(42) 


[ d?x [V • E 2 + + 2/i2$B. ■ (43) 

Thus, ^J-i^Bz is an effective charge which couples to 
the flux of Afi, while Hi^Ba + f^2^Bz is an effective 
charge that couples to the flux of Taking this into 
account, we arrive to the Aharonov-Bohm phase (411. 


In fact, Eqs. (42) and (43) give us a working defini¬ 


tion for the charges of both the Afj, and fields. The 
first equation gives Qa = while the second gives 

0 = —Qz+‘2^^^l^BA+‘^^^2^Bz■ Both of these relations can 
be checked against their definitions given in Sec. IV A 


Notice that according to this definition, Z^ does not 
carry a charge, which is expected since it is a short range 
field and its charge is screened. Now, in order for any 
number of vortices to have a zero net A^^-charge, we must 
have = 0 for the winding numbers violating 

this condition will mean that the system has a logarith¬ 
mically divergent energy. In other words, any collection 
of vortices will be confined if and only if it has a zero net 
charge. Therefore, our vortices are the dynamical realiza¬ 
tion of the Gornalba-Propitius-Wilczek classical confine¬ 
ment phenomenon Haim- 

In this work we did not discuss the stability of our 
vortices as it is beyond the scope of this paper. However, 
we expect the ones with n = ±1 to be stable against de¬ 
cay. As we discussed above, the Abrikosov-Nielsen-Olesen 
and global vortices lie on the opposite sides of the inter¬ 
val fj .1 G [0,oo] (Fig. [^. The stability of both kinds of 
vortices were studied in [52], and it was found that the 
ones with n = 1 are stable, as expected on topological 
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grounds. Indeed, an analysis that follows the lines of [H] 
should be repeated for our vortices to insure their stabil¬ 
ity. However, since n = 1 vortices on the boundaries of 
the Hi interval were found to be stable, it is implausible 
that they lose stability in between as we vary fj-i- 
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